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Abstract 

We achieve a bosonization of one-dimensional ideal gas of exclusion statistics 
A at low temperatures, resulting in a new variant of c = 1 conformal field 
theory with compactified radius R = y/l/X. These ideal excluson gases ex- 
actly reproduce the low-T critical properties of Luttinger liquids, so they can 
be used to characterize the fixed points of the latter. Generalized ideal gases 
with mutual statistics and non-ideal gases with Luttinger-type interactions 
have also similar behavior, controlled by an effective statistics varying in a 
fixed-point line. 
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Recently a new combinatoric rule for many-body state counting [|l],@, which essentially 
is an abstraction and generalization of Yang- Yang's counting ||,|J in Bethe ansatz solvable 
models, is shown to be applicable to elementary excitations in a number of exactly solvable 
models for strongly correlated systems [1,2,4-7] in one, two and higher dimensions. This 
has led to the notion of fractional exclusion statistics [TJ (FES) and associated generalized 
ideal gases [@] (GIG). (For abbreviation, we will call particles or excitations obeying FES as 
exclusons.) An important issue to be addressed is the relevance of FES to realistic gap less 
systems for which state counting is somewhat obscure due to strong correlations among 
particles. In this Letter, we suggest that at least for some strongly correlated systems or 
non- Fermi liquids, their low- energy or low-temperature fixed point may be described by a 
GIG associated with FES, similar to the way that of Landau-Fermi liquids by ideal Fermi 
gas ]|. 

A well-established class of non- Fermi liquids is the Luttinger liquids in ID, proposed by 
Haldane [|9] . To produce in this case a testimony to our above suggestion, we will show that 
the low-T critical properties of Luttinger liquids are exactly reproduced by those of ID ideal 
excluson gases (IEG), if one identifies the controlling parameter |J of the former with the 
statistics A of the latter. So IEG can be used to describe the fixed points of Luttinger liquids. 
A main development here is that we have succeeded in bosonizing ID excluson systems at 



low T, a la Tomonaga | lQfl . It results in a new variant of conformal field theory (CFT) 



with central charge c = 1 and compactified radius |TT| R = Jl/X. The particle-hole duality 



between A and 1/ A in IEG [f||| gives rise to a duality between R and 1/R in this variant of 
c = 1 CFT. 

We have also studied the effects of mutual statistics between different momenta and 
Luttinger-type (density-density) interactions among exclusons. In both cases, the low-T 
behavior is controlled by an effective statistics X e ff for excitations near the Fermi points, 
the same way that of IEG by A. In ID both the momentum independent part of interactions 
and change in chemical potential [i are relevant perturbations P,|T2"|], resulting in a continuous 
shift in the fixed-point line parametrized by A. 



Consider a GIG of iVo particles on a ring with size L. Single-particle states are labeled 
by hi. The total energy and momentum are given by E = kf and P = J^ki- We assume 
that in the thermodynamic limit the density, p a {k), of available single-particle states is 
linearly dependent on the particle density, p(k). By definition, the statistics matrix is given 
by the derivative 

g(k,k') = -5p a (k)/5p(k>). (1) 

The system is called an IEG of statistics A (with no mutual statistics between different 
momenta), if g(k,k') = X5(k — k'), or p a {k) + Xp{k) = po(k) [|J, where po(k) = l/2n is the 
bare density of single-particle states. Thus, A = 1 corresponds to fermions, and A = to 
bosons. The thermodynamics of IEG is shown |2| to be determined by the thermodynamic 
potential 

Q = / dk ln(l + w(k, T) _1 ), (2) 

2% J-oo 

with the function w{k,T) = p a (k)/p(k) satisfying 

w(k, T) A [1 + w(k, I 1 )] 1 "* = e {k2 ~^ /T . (3) 

In the ground state, there is a Fermi surface such that p(k) = l/lirX for \k\ < /cj? and p(k) = 
for | A; | < kp. Then the Fermi momentum is given by kp = nXdo, and the ground state energy 
and momentum by E /L = 7t 2 A 2 Jq/3, P — 0, with the average density d = N /L. 

Now let us examine possible excitations in the IEG. First there are density fluctua- 
tions due to particle-hole excitations, i.e. sound waves with velocity v s = vp = 2kp (see 
below). Moreover, by adding extra M particles to the ground state, one can create particle 
excitations, and by Galileo boost a persistent current. Our observation is that the veloc- 
ities of these three classes of elementary excitations in IEG satisfy a fundamental relation 
that Haldane years ago used to characterize the Luttinger liquid 0. Indeed, shifting iVo to 
N = N + M, the change in the ground state energy is SiE = ir(Xkp)M 2 , while a persistent 
current, created by the boost of the Fermi sea k — > k + ttJ/L, leads to the energy shift 
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82E0 = itikp/X) J 2 - (Due to periodic boundary conditions, M and J are constrained by 
M = J (mod 2).) Therefore the total change in energy and in momentum, due to charge 
and current excitations are 

5E = ^(v N M 2 +vjJ 2 ), 5P = 7T(d + ^-)J, (4) 

respectively, with 

V N = V S X, Vj = V s /X, V s = y/V N Vj. (5) 

These coincide with the well-known relations |J in Luttinger liquid theory, if we identify 
A with Haldane's controlling parameter exp(— 2(f). So it is interesting to see whether all 
critical exponents of Luttinger liquids would be reproduced simply by IEG. 

To calculate the exponents in IEG, we need to develop a bosonization approach. Follow- 
ing Yang and Yang fljlll, we introduce the dressed energy e(k,T) by writing 



w(k,T) = e« k ' T ^ T . (6) 

The point is that the grand partition function Zq, corresponding to the thermodynamic 
potential (2), is of the form of that for an ideal system of fermions with a complicated, In- 
dependent energy dispersion given by the dressed energy: Zq = Y\k{^-+^~ t<yk,T ^ T ). However, 
this fermion representation is not very useful, because of the implicit T-dependence of the 
dressed energy. We observe, nevertheless, that in the low-T limit, the T-dependence of 
e(k,T) can be ignored: according to eqs. (3) and (6), e(k,T) = e(k) + 0(e _ ' e '/ T ), where 

, M ( (k 2 - kj 7 )/X + kjr, \k\ < k F , 

e(k) = (7) 
I k 2 , \k\ > k F . 

Thus, the low-T grand partition function can be obtained from the effective Hamiltonian 
given by i? e // = J2k e (k) c\ck, where c\ are fermionic creation operators. 

Another simplification in the low-T limit is that we need to consider only low-energy 
excitations near Fermi points k ~ ±/c.f, where the left- and right-moving sectors are separable 
and decoupled: H e ff = H + + H_, and H± has a linearized energy dispersion 



f ±v F {k =f k F ) + kp, \k\>k F , 
e±{k) = i (8) 

^ ±v F (k =f + kp, \k\ < k F . 

We note the 'refractions' at k = ±k F . In spite of this peculiarity, we have succeeded in 
bosonizing the effective Hamiltonian as follows. The density fluctuation operator at k ~ k F 
is constructed as follows 

P ( q +) = H : Ck-qCk ■ + '■ C k+Xq C k ■ 

k>kp k<kp—\q ^g-j 

+ V] '■ (^k-k F Cfc : 

— T^^+fcFH-? 

kp-Xq<k<kp A 

for g > 0. A similar density operator p^"- 1 can also be defined for k ~ — kp. Within the 
Tomonaga approximation ||14|| , in which commutators are taken to be their ground-state 
expectation value, we obtain 

[P q ± \p { t ) '} = ^K^ [H ± , pW) = ±v F qp q . (10) 
which describe ID free phonons with the sound velocity v s = v F . Introducing normal- 



ized boson annihilation operators b q = y2n/ XqL p q + \ b q = yz7T /XqLp q ^, the bosonized 
Hamiltonian satisfying (|10"D is given by 

He = v s {Y,q(b% + b\b q ) + \j[XM* + (11) 

The construction of the bosonized momentum operator is a bit more tricky. Each term 
in eq. (9) should carry same momentum q, therefore the fermion created by c k carries a 
dressed momentum p, which is related to k by 

k — k F + (k F /X), k > kp, 
p(*0 = { k/X, \k\<k F ,. (12) 

k + kp — (kp/X), k < —k F . 

In terms of this variable, the linearized dressed energy e(p) is of a simple form: e±(p) = 
±v s (p =F p F ) + p, with p F = kp/X. The bosonized total momentum operator, corresponding 
to the fermionized P = J2kP{k) c \ c k, is 



P = E - ~ b %) + <d + M/L) J. (13) 

q>0 



In the coordinate-space formulation, the normalized density field p(x) is given by p(x) 
Pr{x) + Pl(x): 

M 



Pr(x) = Tpr + E ^q/2*L\(e^b q + e^ x b\), (14) 

and Pl{x) is similarly constructed from b q and K. The boson field (fr(x), which is conjugated 
to p(x) and satisfies [<f>(x), p(x')] = i5(x — x'), is 0(x) = 4>r(x) + 0l(x) with 

and a similar 4>l{x). Here M and J are operators with integer eigenvalues, and 0o is an 



angular variable conjugated to M: [0 O , M] = i. The Hamiltonian (|IT ) becomes 



# B = ^ [n(x) 2 + (<9 x .X(x)) 2 ], (15) 

Z7T JO 

where n(x) = vrA 1 / 2 p(a;) and X(x) = A~ 1 / 2 0(x). With X(x,t) = e iHt X(x)e~ iHt , the La- 
grangian density reads 

C = ^d a X(x,t)d a X(x,t). (16) 

Z7T 



We recognize that C is the Lagrangian of a c = 1 CFT [|TTJ. Since 0o is an angular 
variable, there is a hidden invariance in the theory under <fi — > + 27r. The field X is thus 
said to be "compactified" on a circle, with a radius that is determined by the exclusion 
statistics: 

X ~ X + 2^, i? 2 = 1/A. (17) 

States V[X]|0) or operators V[X] are allowed only if they respect this invariance, so quantum 
numbers of quasiparticles are strongly constrained. The particle-hole duality [f|||, i.e. A <-» 
1/A and M J, in eq. fly]) is just a duality R <-» 1/R. Moreover, the partition function of 
IEG (in the low-T limit) can be rewritten as Z = Tr^[q L o q L °], where q = e Ws ^ T . The zeroth 
generators of the Virasoro algebra are Lq' l = v' 1 !!^ + fi[JR =F M/R] 2 . The constraint 



M = J (mod 2) makes our variant of c = 1 CFT have an unusual spectrum and duality 
relation fllTf . 

We also see that the Hamiltonian flI5|) agrees with Haldane's harmonic fluid description 
of Luttinger liquids [|15j. Thus the critical properties of IEG reproduce those of Luttinger 
liquids. Or one may say that IEG can be used to characterize the fixed points of Luttinger 
liquids. 

Using bosonization techniques, the computation of correlation functions and critical ex- 
ponents in IEG is similar to that for Luttinger liquids |]9|JT5|. A careful construction for the 
allowed boson field with charge- 1 leads to 



oo 



^ B (x,t) = p(x)^ 2 £ e i0 ™ : e^WA 1 / 2 )^-) . 

m=— oo 

. e i(X 1 / 2 -2m/X 1 / 2 )X L (x+) . ^ ( lg ) 

where the hermitian, constant- valued operators O m satisfy [O m ,O m i} = m(m — m') ||16|| . 
The excluson operator reads ^\(x) =: ty B (x)e % * x 2 ( Xr ^~ Xl ^ :, obeying ^\(x)^\(x') — 
e «7rAsgn(x-a;')^>t ( x ')ty x (x) = for x x' . The multi-sector density operator for exclusons is 

p(x) = ^{(x)^x(x) = ¥ B (x)^ B (x) 
= p(x)Y, m ■ exp{t2m[X R (x) - X L {x)}/ X 1 ' 2 } : (19) 
Dynamical correlation functions can be easily calculated: 

(p(x,t)p(o,o)> « 4 

oo 

+ J2 A m i^jpTA cos(2nd mx) 

m=l 

G(x,t;X) = (*t(x,t)* A (0,0)) 

oo 

m=— oo - + 

with A m and B m regularization-dependent constants. These correlation functions coincide 
with the asymptotic ones in the Calogero- Sutherland model. 

The single-hole state, i.e. \l/jy A |0) = ^\(X — ► 1/A)|0), with charge — 1/A alone is not 
allowed. The minimum allowed multi-hole state is given by ^f\, x (xi)...^\/ x (x p )\0) if A = p/q 
is rational. One may obtain, e.g. , 



1 / 1 J_ 

+ (2vrJ ) 2 A [x 2 + + x 2 _ 



([*I/a(z, *)F[*i/a(0, 0)H ~ [G(z, t; 1/A)f. (20) 

A more interesting allowed operator is what creates q particle excitations accompanied by 
p hole excitations: n(x,t) = [9^(x,t)] q [W 1 , x (x,t)] p . We note the similarity of this operator 
to Read's order parameter |17j for fractional quantum Hall fluids (in bulk). Its correlation 
function can be calculated by using Wick's theorem: 

(h(x, t)n(0, 0)) ~ [G(x, t; \)} q [G{x, t; 1/A)f. (21) 

If the contribution from the m = sector dominates, then one gets (n(x, t)n(0, 0)) ~ 
(x - v s t)- {p+q \ 

Now we turn to discussing the effects of mutual statistics. Consider a GIG with the 
statistical matrix ([[]) given by g(k — k') = 5(k — k') + $(/c — k'). Here = $(— k) is a 
smooth function and stands for mutual statistics between particles with different momenta, 
in contrast to IEG, for which <$>(k) is a 5-function. The thermodynamic properties of GIG 
is also given by eq. (J2|), but now w(k,T) satisfies instead an integral equation 0,f|] which, 
in terms of the dressed energy (^), is of the form 

e(k, T) = e (k) + T J ^ S(Jfe - k') ln(l + e"^' T ^ T ) , 

where e (k) = k 2 , and we have shifted the dressed energy by chemical potential //. At T = 0, 
the Fermi momentum kp is determined by e(±fcp) = 0. Introduce 

— a(k)P(k), (22) 

^— $(k - k') a(k') . (23) 
-k F 2n 

Then both p{k) and e(fc) in the ground state satisfy an integral equation like 

a(k) — a (k) — ($ • a) (k; —k F , k F ) • (24) 

The dressed momentum p(k) is related to p(k) by p'(k) = 2irp(k) and p(k) = —p(—k). 
The ground state energy is given by E /L = (eo • p)[— kp, kp] = (e • po)[— kp,kp}. These 
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equations are of the same form as those in the thermodynamic Bethe ansatz J3], hence the 
Luttinger-liquid relation v s = y/vNVj remains true ]TjJ. A simple proof can be sketched 
as follows. The sound velocity is well-known: v s = de(pp) /dpp. The charge velocity is 
given by vn = v s z(kp)~ 2 , where the dressed charge z(k) is given by the solution to the 
integral equation z(k) = 1 — ($ • z)(k; —kp,kp}. This relation can be easily derived from 
the definitions vn = Ldp/dN and z{k) = —5e(k)/5p. To create a persistent current, let us 
boost the Fermi sea by ±kp — > ±kp + A, where A = z(kp) j 'Lp{k,p). Then the total energy 
of the state with the persistent current is 

E A /L = (e • p A ) [-k F + A,k F + A] 

= (e A -p )[-k F + A,k F + A], (25) 

where p A (k) = p (k)-($-p A )(k; -k F +A, k F +A] and e A (k) = e Q (k)-($-e A )(k] -k F +A, k F + 
A]. Now, using the last expression of E A , it is easy to show E A — E = LA 2 e'(k F )p(k F ) = 
(27r / ' L)v s z(k F )~ 2 . This verifies v j = v s z{k F ) 2 . In view of eq. (H), at low energies GIG looks 
like IEG with 

A e// = z{k F )~ 2 . (26) 
It can be shown that it is the effective statistics fl2"E|) that controls the low-T critical 



properties of GIG, as A does for IEG. Linearization near the Fermi points and bosonization 
of the low-energy effective Hamiltonian go pretty the same way as before for IEG. The 
only difference now is that the linearized dispersion for dressed energy e±(k) = ±e'(k F )(k =F 
k F ) + p = ±v s (p(k) TPf) + is smooth at k ~ ±k F . So bosonization is standard and the 
bosonized Hamiltonian is the same as eq. ( PD for IEG, only with A replaced by A e //. An 
(allowed) ^\ describes the particle excitations near the Fermi surface with both anyon 
and exlcusion statistics being X e ff- in this sense, one may say that the effect of mutual 
statistics is to renormalize the statistics. 

Here we remark that in IEG, $(A;, k') = (A — 1) 5(k — k') is not smooth, so the dressed 
charge has a jump at kp\ z(kp) = 1 and z(kp) = A -1 for kp = kp ± + . The general 
Luttinger-liquid relation is of the form 
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v N = v s [z(kp)z(k F )] \ vj = v s z(kp)z(k F ). (27) 

Finally we examine non-ideal gases, e.g., with general Luttinger-type density- density 
interactions 

H i = T YX U *(PiP\ + PiP\) + VMpI + ~P q P% (28) 
<?>o 

which is added to the Hamiltonian describing a GIG. After bosonization, the total Hamil- 
tonian remains bilinear in densities, so it is trivial to diagonalize it using the Bogoliubov 
transformation. The diagonalized Hamiltonian is again of the harmonic-fluid form ([LTD 
with b q and b q replaced by corresponding operators for Bogoliubov quasiparticles, and with 
the velocities renormalized: v s — > v s = \(v s + U ) 2 — Vq 2 ) 1 / 2 , vn — > = v s e~ 2 ^ , and 
v j ~^ Vj = v s e 2l ^°. Thus, the Luttinger-liquid relation (|5|) survives, with A e // of GIG renor- 
malized to 

A e// = e- 2 *\ with tanh(2^ ) = - ~ ~ 7 • ( 29 ) 

Note that the new fixed point depends both on the position of the Fermi points and on the 
interaction parameters U$ and Vo, leading to "non-universal" exponents. 

In passing we observe several additional implications of this work: 1) Our bosonization 
and operator derivation of CFT at low energies or in low-T limit can be applied to Bethe 
ansatz solvable models, including the long-range (e.g. Calogero-Sutherland) one. 2) Here 
we have only consider one-species cases, i.e. with excitations having no internal quantum 
numbers such as spin. Our bosonization and characterization of Luttinger liquids are gener- 
alizable to GIG with multi-species, presumably with the effective statistics matrix related to 
the dressed charge matrix. 3) The chiral current algebra in eq. ( fLOD with A = 1/m coincides 



with that derived by Wen |19[ for edge states in v = 1/m fractional quantum Hall fluids. 
So these edge states and their chiral Luttinger-liquid fixed points can be described in terms 
of chiral IEG. 

In conclusion, we have shown that ID IEG without mutual statistics can exactly repro- 
duce the low-energy and low-T properties of (one-component) Luttinger liquids. Moreover, 
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mutual statistics and Luttinger-type interactions in a GIG only shift the value of A e //. Thus 
the essence of Luttinger liquids is to have an IEG obeying FES as their fixed point. It 
is conceivable that some strongly correlated systems, exhibiting non-Fermi liquid behavior, 
in two or higher dimensions can also be characterized as having a GIG with appropriate 
statistics matrix as their low-energy or low-temperature fixed point. 
This work was supported in part by NSF grand PHY-93094598. 



11 



REFERENCES 

[1] F. D. M. Haldane, Phys. Rev. Lett. 67, 937 (1991); and in Proc. 16th Taniguchi Sym- 
posium, eds. N. Kawakami and A. Okiji, Springer Verlag (1994). 



Y. S. Wu, Phys. Rev. Lett. 73, 922 (1994). 

C. N. Yang and C. P. Yang, J. Math. Phys. 10, 1115 (1969). 

D. Bernard and Y. S. Wu, in Proc. 6th Nankai Workshop, eds. M.L. Ge and Y.S. Wu, 
World Scientific (1995). 

C. Nayak and F. Wilczek, Phys. Rev. Lett. 73, 2740 (1994). 
Z.N.C. Ha, Nucl. Phys. 435, 604 (1995). 



[2] 
[3] 

K 

[5] 
[6] 

[T 
[8 
[9 
[10 

[11 

[12 

[13] B. Sutherland, J. Math. Phys. 12, (1971) 251. 

[14] Alternatively, we may extend the Fermi sea to that for two separate right- and left- 
moving fermions with linear dispersion for k covering the whole real axis in each sector. 
Then our commutation relations are exact. (See D.C Mattis and E.H. Lieb, J. Math. 
Phys. 6, 304 (1965).) 

[15] F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 



Y. Hatsugai, M. Kohmoto, T. Koma and Y.S. Wu, preprint |cond-mat/ 9504 112 
See, e.g. , R. Shankar, Rev. Mod. Phys. 66, 129, (1994). 
F. D. M. Haldane, J. Phys. C 14, 2585 (1981). 
S. Tomonaga, Prog. Theor. Phys. 5, 544 (1950). 

See, e.g. , R. Dijkgraaf, E. Verlinde and H. Verlinde, Commun. Math. Phys. 115, 649 
(1988). The standard duality in usual c = 1 CFT is R <-> 2/R. 



See, e.g., H. J. Schulz, preprint |cond-mat/9503150 . 



12 



[16] Here operators O m give rise to the Klein factor necessary for correct commutation 
relation for and ^ B . 

[17] N. Read, Phys. Rev. Lett., 62, 86 (1988). 

[18] F. D. M. Haldane, Phys. Lett. 81 A, 153 (1981); N. M. Bogolyubov, A. G. Izergin and 
V. E. Korepin, J. Phys. A 20, 5361 (1987). 

[19] X.G. Wen, Phys. Rev. Lett. 64 , 2206 (1990). 



13 



